hep-ph/ 0307091 
SHEP/0320 



O 

o 

(N 



(N 
> 



a* 



Lepton Flavour Violation from Yukawa Operators, 
Supergravity and the See-Saw Mechanism 

S. F. King and I. N. R. Peddie 



Department of Physics and Astronomy, University of Southampton 

(N 

<N ■ Southampton, SO 17 IB J, U.K 



Abstract 



We investigate the phenomenological impact of different sources of lepton flavour vi- 
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O ' olation arising from realistic models based on supergravity mediated supersymmetry 

o breaking with Yukawa operators. We discuss four distinct sources of lepton flavour 



violation in such models: minimum flavour violation, arising from neutrino masses 
and the see-saw mechanism with renormalisation group (RG) running; supergravity 



^ ■ flavour violation due to the non-universal structure of the supergravity model; flavour 

5-H ' 



violation due to Froggatt-Nielsen (FN) fields appearing in Yukawa operators devel- 
oping supersymmetry breaking F-terms and contributing in an non-universal way to 
soft trilinear terms; and finally heavy Higgs flavour violation arising from the heavy 
Higgs fields used to break the unified gauge symmetry which also appear in Yukawa 
operators and behave analagously to the FN fields. In order to quantify the relative 
effects, we study a particular type I string inspired model based on a supersymmetric 
Pati-Salam model arising from intersecting D-branes, supplemented by a U(l) family 
symmetry. 
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1 Introduction 



Lepton flavour violation (LFV) has been long known to be a sensitive probe of new 
physics in supersymmetric (SUSY) models [UE]- LFV arises in SUSY models due to 
off-diagonal slepton masses in the super-CKM basis in which the Yukawa matrices are 
diagonal. Such flavour violation could arise either directly at the high energy scale 
due to primordial string or SUGRA effects, or be generated radiatively by the renor- 
malisation group equations, for example in running a grand unified theory (GUT) 
from the Planck scale to the GUT scale (due to the presence of Higgs triplets) |Hj or 
in running the minimal supersymmetric standard model (MSSM) with right-handed 
neutrinos from the Planck scale down to low energies, through the scales at which 
the right-handed neutrinos decouple. 

Even in minimal supergravity (mSUGRA), where there is no flavour violation at 
the high energy scale, the presence of heavy right-handed neutrinos as required by the 
see-saw mechanism explanation of small neutrino masses will lead inevitably to LFV 
jH Ej. The recent neutrino experiments which confirm the matter enhanced Large 
Mixing Angle (LMA) solution to the solar neutrino problem [0], together with the 
atmospheric data [Zj, show that neutrino masses are inevitable, and, assuming SUSY 
and the see-saw mechanism, hence show that LFV must be present. For example this 
has recently been studied in mSUGRA models with a natural neutrino mass hierarchy 
jH]. There is in fact a large literature on this subject [§]. 

Recently it has been realised that in realistic string inspired models based on su- 
pergravity mediated supersymmetry breaking, in which the origin of Yukawa matrices 
Yij is due to Froggatt-Nielsen (FN) operators fU] of the form ~ 8 niJ , where is 
an integer power, there may be a new and dangerous source of LFV which arises when 
the FN fields 9 develop supersymmetry breaking F-terms F$ ~ ^3/2$ leading to non- 
universal soft trilinear terms AAij = FgdglnYij which implies AAy ~ nijm 3 / 2 
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[T^l IT3] . The effect is independent of the vacuum expectation value (vev) of the FN 
field 9, and is present even in minimum flavour violation scenarios such as mSUGRA. 

In this paper we shall explore the phenomenological impact of the new source of 
LFV arising from FN fields discussed above, and compare it to the more usual sources 
of LFV arising from right-handed neutrinos, and non-universal SUGRA models in 
order to gauge its relative importance. A phenomenological analysis is necessarily 
model dependent, and so we shall study a particular type I string inspired model based 
on a supersymmetric (SUSY) Pati-Salam model arising from intersecting D-branes, 
which was introduced in However in order to explore the effects of interest, it 
is necessary to supplement this model by a U(l) family symmetry, and introduce FN 
fields so as to provide a realistic description of quark and lepton masses and mixing 
angles, including those of the neutrino sector. Recently a global \ 2 analysis of a 
realistic SUSY Pati-Salam model was performed ^H], arid a good fit to the quark 
and lepton mass spectrum was obtained based on a FN operator analysis with a U(l) 
family symmetry. It is therefore natural to combine the models in and ^Hj in 
order to provide a realistic framework for studying the new LFV effects arising from 
the FN fields, and to compare this to the effects on non-universal SUGRA and also 
right-handed neutrinos in a model that gives a good fit to the neutrino data. 

Of course in combining the two models we are taking some liberties with string 
theory. In particular we assume that the combined model corresponds to the low 
energy limit of a string model as in [2] , but with the addition of an extra state, which 
is a Froggatt-Nielsen ^U] family field, 9. We assume that since the model without 9 
can be extracted from a string model, then so can the model with 9, but we make no 
attempt to derive it. We emphasise that the main motivation for combining the two 
approaches is to explore the phenomenology of LFV in a "realistic" framework. One 
by-product of doing this is that we identify a genuinely new source of LFV that has 
not been considered at all in the literature, namely the heavy Higgs fields H that break 
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the unified gauge symmetry at high energies. These heavy Higgs fields also appear 
in the operators which describe the Yukawa couplings, and they can be expected 
to behave in a similar way to the FN fields 9, and give rise to LFV analagously. 
The combined model has a number of attractive features: it includes approximate 
third family Yukawa unification, the number of free parameters is restricted to eight 
undetermined free parameters related to supergravity, and the model gives a good fit 
to all quark and lepton masses and mixing angles. 

In order to study the phenomenological effect of the different sources of LFV, we 
generalise the Goldstino Angle parametrisation of the dilaton and moduli fields S, Tj 
to include a parametrisation of the SUSY breaking F-terms for the FN fields 9 and 
heavy Higgs fields H. There are four distinct sources of lepton flavour violation in 
this model: minimum flavour violation, arising from neutrino masses and the see-saw 
mechanism with renormalisation group (RG) running; supergravity flavour violation 
due to the non-universal structure of the supergravity model; FN flavour violation 
due to the FN fields developing supersymmetry breaking F-terms and contributing in 
an non-universal way to soft trilinear terms; and finally heavy Higgs flavour violation 
arising from the Higgs fields used to break the unified gauge symmetry which may 
behave analagously to the FN fields. We propose four benchmark points at which 
each of these four sources separately dominate. We then perform a detailed numerical 
analysis of LFV arising from the four benchmark points. The numerical results show 
that LFV due to FN fields is the most sensitive source in the sense of leading to larger 
limits of m 3 / 2 , however we find that the gluino mass is relatively light in these cases 
which tends to reduce fine-tuning. We also find that in some cases the LFV effects 
from Yukawa operators in the presence of the seesaw mechanism can be less than 
without the seesaw mechanism. 

The outline of this paper is as follows. In section El we introduce the specific 
model that we shall study, discuss the symmetries of the model, and the Yukawa 
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and Majorana operators, and for particular choices of the order unity coefficients, 
show that this leads to a good fit to the neutrino data, with a prediction for the 
unmeasured 6* 13 . In section|3]we discuss the soft SUSY breaking aspects of the model. 
We parametrise the SUSY breaking F-terms, give the soft scalar masses, including 
the D-term contributions, give the soft gaugino masses and soft tirlinear masses, and 
explain why these are expected to lead to large flavour violation. In section 0] we 
give the results of a numerical analysis of the model, focussing on four benchmark 
points designed to highlight the four different sources of LFV. Finally we present our 
concluding remarks in section |S] 

2 The Model 

2.1 Symmetries and Symmetry Breaking 

The model defined in Table ^ is an extention of the Supersymmetric Pati-Salam 
model discussed in ref . [T2] , based on two D5 branes which intersect at 90 degrees and 
preserve SUSY down to the TeV energy scale. The string scale is taken to be equal 
to the GUT scale, about 3 x 10 16 GeV. 

The extension is to include an additional U(1)f family symmetry and the FN 
operators as in ^H] (see also (El)- The present 42241 Model is then just the 4224 
Model of Jl] augmented by a U(l)j? family symmetry. The purpose of this extension 
is to allow a more realistic texture in the Yukawa trilinears Y a b c , along the lines of the 
recent operator analysis in [To] . 

The quark and lepton fields are contained in the representaions F, F which are 
assigned charges Xp under XJ(l)p. In Table list two equivalent sets of charges 
U(1)f and U(l)^, where \J(1)f is anomaly free, but U(l)-^ is equivalent for all prac- 
tical purposes and has much simpler charge assignments. The field h represents both 
Electroweak Higgs doublets that we are familiar with from the MSSM. The fields H 
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Table 1: The particle content of the 42241 model, and the brane assignments of the 
corresponding string 



and H are the Pati-Salam Higgs scalars; 1 the bar on the second is used to note that 
it is in the conjugate representation compared to the unbarred field. 

The extra Abelian U(l)p gauge group is a family symmetry, and is broken at the 
high energy scale by the vevs of the FN fields fU] 0,9, which have charges —1 and +1 
under U(l)p, respectively We assume that the singlet fields 6, 9 arise as intersection 
states between the two D5— branes, transforming under the remnant U(l)s in the 
4224 gauge structure. In general they are expected to have non-zero F-term vevs. 

The two SU (4) gauge groups are broken to their diagonal subgroup at a high scale 
due to the assumed vevs of the fields ipi, y?2 |14j . The symmetry breaking at the scale 
M x 

SU(4) ® SU(2) L ® SU(2) H -> SU(3) <g> SU(2) L ® U(l)y (1) 

1 We will also refer to these as "Heavy Higgs"; this has nothing to do with the MSSM heavy 
neutral higgs state H° 
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is achieved by the heavy Higgs fields H, H which are assumed to gain vevs [T5j ) 

(H ab ) = (u H ) = V5 a J b 2 ~ M x ; (H ax ) = (V H ) = V5 A J 2 X ~ M x (2) 

This symmetry breaking splits the Higgs field h into two Higgs doublets, h\, h 2 . Their 
neutral components then gain weak-scale vevs 

(h° 1 ) = v 1 ; (h° 2 ) =v 2 ; tfw/3 = v 2 /v 1 . (3) 

The low energy limit of this model contains the MSSM with right-handed neutrinos. 
We will return to the right handed neutrinos when we consider operators including 
the heavy Higgs fields H, H which lead to effective Yukawa contributions and effective 
Majorana mass matrices when the heavy Higgs fields gain vevs. 

2.2 Yukawa Operators 

The Yukawa operators, responsible for generating effective Yukawa couplings, have 
the following structure 2 |16j : 

where the integer p(i,j) is the total U(1)f charge of Fi + Fj + h and HH has &U(1)f 
charge of zero. The tensor structure of the operators in Eq0]is 

= F aa F Px hlH lz HP w e p ^ (5) 

One constructs [in] SU(4)p5 invariant tensors C^J that combine 4 and 4 representa- 
tions of SU(4)p5 into 1, 6, 10, 10 and 15 representations. Similarly we construct 
SU(2)^ tensors B^ z w that combine 2 representations of SU(2) into singlet and triplet 
representations. These tensors are contracted together and into to create 

singlets of SU(4)ps, SU(2) L and SU(2)#. Depending on which operators are used, 

different Clebsch-Gordan coefficients (CGCs) will emerge. 

2 We note that due to the allocation of charges, and since the effective Yukawa operators include 
the fields FjFjh with overall charge positive, the field 6 does not enter the Yukawa operators. 
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We look at two different models for the Yukawa sector, which we refer to as model I 
and model II. The models represent different 0(1) paramaters a, a', a" in the following 
operator texture [To] : 



O 



{a n O Fc + a![ x O" Ae )£ (a 12 Ee + a' 12 0' Cb )e 3 (a' 13 ,c/ + a'{ 3 0" Ee )e 



(a 21 Dc )e 4 
(a 3 iO Fc )e 4 



{a 22 Bc + a' 22 O lF f)e 2 (a 23 Ee + a' 23 Q 
(a 32 Ac + a'O lFe )e 2 a 33 



/Bc\ 



(6) 



where the operator nomenclature is defined in Appendix O For convenience, from 
this point on, we define: 

S 



EH 



M\ 



(7) 



and 
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M 



x 



EqElthen yields the effective Yukawa matrices 
Y U (M X ) = 








22 5^° 6 



n' 2A2.2 

a 3 2 5 o e 



^33 



Y d {M x ) 



\ 2 \plb 



2^3 



«2i^ 4 (a 22 JlS + a' 22 ^)e* 



Y e (M_ 



x, 



a 3 i|5e 4 
-a 31 Ue 4 



a 32 



V25e 2 



a' -i. 

A33 



-2 2 3vf4 + a' 22 ^)^ 



a 23 v / 25e 2 



2X2 



- a 23VF 
%3 



Y V [M_ 



x, 



» 2 V26 





3^5 



a 12 25e 3 



a» 3 5 3 e 



A 22 5 y5 

a' 

u 32 5° e 



5V a 23 2S 



«33 



(9) 



(10) 



11] 



(12) 



The order unity coefficients a^, a'^ of the operators are adjusted to give a good 
fit to the quark and lepton masses and mixing angles, and take the values given in 
Table El Note that the two models differ only in the choice of a" 3 , which is taken to be 
zero in model I. Model I consequently has a lower rate for \x — > e^y, and model II has 
a higher \x — > rate due to the non-zero 13 element of the neutrino Yukawa matrix, 
as can be understood from the analytic results in jH]. The fits assume 5 = e = 0.22. 
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Model II 
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0.55 


an 


-0.92 


-0.92 


au 


0.33 


0.33 
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1.67 


1.67 


^22 


1.12 


1.12 


^23 


0.89 


0.89 


031 


-0.21 


-0.21 
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2.08 


2.08 


Ol2 


0.77 


0.77 


0'l3 


0.53 


0.53 


°22 


0.66 


0.66 


°23 


0.40 


0.40 


«32 


1.80 


1.80 


<X 


0.278 


0.278 


<3 


0.000 


1.000 


An 


0.94 


0.94 


A X2 


0.48 


0.48 


A13 


2.10 


2.10 


^22 


0.52 


0.52 


^23 


1.29 


1.79 


^33 


1.88 


1.88 



Table 2: The a, a' and a" paramaters for model I and model II 
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2.3 Majorana Operators 



We are interested in Majorana fermions because they can contribute neutrino masses 
of the correct order of magnitude via the see-saw effect. The operators for Majorana 
fermions are of the form 

There do not exist renormalisable elements of this infinite series of operators, so 
n < 1 Majorana operators are not defined 3 . A similar analysis goes through as for 
the Dirac fermions; however the structures only ever give masses to the neutrinos, 
not to the electrons or to the quarks. 4 

It should be noted that the Majorana neutrinos will not affect the A-terms, as 
these operators do not contribute to the Yukawas. The RH Majorana neutrino mass 
matrix is: 

M RR (M X ) 



M- 



:>,:', 



A u Se 8 A 12 Se e A l3 5e A 



A 12 5e e A 22 5e* A 23 5e 2 
A 13 5e* A 23 5e 2 A 33 



(14) 



2.4 Neutrino sector results 

The neutrino Yukawa matrix in EqJ12land the heavy Majorana mass matrix in EqfTH 
imply that the see-saw mechanism satisfies the condition of sequential dominance |17j . 
leading to a natural neutrino mass hierarchy mi <C m 2 <C m 3 with no fine-tuning. 
The dominant contribution to the atmospheric neutrino mass m 3 comes from the 
third (heaviest) right-handed neutrino, with the leading subdominant contribution to 
the solar neutrino mass m 2 coming from the second right-handed neutrino. In such 



3 Except for the 33 neutrino mass term; this is allowed because of string theoretic effects 
4 To see this note that the form of the two H vevs is symmetric, and proportional to 626% ■ 
Symmetric structures will then contract to give neutrino mass terms. Antisymmetric structures will 
contract to give zero. As any structure can be written as a sum of a symmetric and an antisymmetric 
part, we see immediately that the only mass terms can be given to the neutrinos because of the form 
of the vevs in Eq. (J2J 
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a natural scenario, the large atmospheric angle is due to the large ratio of dominant 
neutrino Yukawa couplings tan 623 ~ ^23/^33; an d the large solar angle is due to the 
large ratio of leading subdominant Yukawa couplings tan 612 ~ \[2Y^ 2 j '(Y 22 — Y^ 2 ). 

Model I for the Yukawa sector is taken from a global analysis of a SUSY Pati- 
Salam model enhanced with an Abelian flavour symmetry [T^]. At one-loop order the 
Yukawa runings only depend on the other Yukawas and the gauge couplings. Since 
Model II only differs from Model I in the neutrino Yukawa, we do not expect the 
quark masses or mixing angles to be different. We also do not expect the charged 
lepton masses to differ by much. 



Observable 


Model I 
Prediction 


Model II 
Prediction 


Experimental 
Values 


sin 2 612 
sin 2 623 
sin 2 9 n 


0.316 
0.553 
5.18 • 10~ 3 


0.308 
0.552 
5.20 ■ 10~ 3 


0.28 ±0.05 
0.50 ±0.15 
< 0.03 


Am ltm 


1.32 ■ 10~ 3 
6.05- 10" 5 


1.33- 10~ 3 
5.91 • 10~ 5 


(2.5±0.8)10~ 3 
(7.0±3.0)10" 5 



Table 3: The neutrino mass differences and mixing angles in model I, model II and 
the experimental limits 

The possibility remains open, however, that the new operator in the 13 Yukawa 
elements could predict either a mass difference or a neutrino mixing angle in violation 
of the results from the various neturino experiments [HI 13 EH QS|- As such, we checked 
our predictions for the mass-differences and the mixing angles for both models, in 
comparison to experiment. The results of this are summarised in Table El 

We note that in both model I and model II, we are within the constraints on 
the neutrino sector. In fact Model II is slightly closer to the central values of three 
observable paramaters (the solar and atmospehric neutrino mixing angles, and the 
atmospheric mass difference). In both cases we predict values of #13 below the current 
limit. 
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3 Soft Super symmetry Breaking Masses 

3.1 Supersymmetry Breaking F-terms 

In jTJj it was assumed that the Yukawas were field-independent, and hence the only 
F-vevs of importance were that of the dilaton (S), and the untwisted moduli (T l ). 
Here we set out the paramaterisation for the F-term vevs, including the contributions 
from the FN field 9 and the heavy Higgs fields H, H. Note that the field dependent 
part follows from the assumption that the family symmetry field, 9 is an intersection 
state. 

F s = V3m 3/2 (S + S)x s (15) 
F Ti = V^ms/a^ + Ti)^ (16) 
Fjjab = s/3m 3/2 H ab (S + Sy X H (17) 
Fjj ax = V3m 3/2 H ax (T 3 + T 3 f Xjj (18) 
F e = VSm 3/2 9(S + sf(T 3 + T 3 fx e (19) 

We introduce a shorthand notation: 

F H H = F Hab H ab ; FjjH = ^F-jj ax H ax . (20) 

ab olx 

3.2 Soft Scalar Masses 

There are two contributions to scalar mass squared matrices, coming from SUGRA 
and from D-terms. In this subsection we calculate the SUGRA predictions for the 
matricies at the GUT scale, and in the next subsection we add on the D-term contri- 
butions. 

The SUGRA contributions to soft masses are detailed in Appendix El From 
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Eq. we can get the family independent form for all scalars: 



m, 



3/2 



(21) 



where 



m 2 R 



m. 



in 
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in 



H 



a 
b L 
b R 



in 



in 



in 



in 



3/2 
2 

3/2 
2 

3/2 
2 

3/2 



>li 



;i - 3^1) 



3J^ 



\a2 ^ 



1 - 3X| 2 



(22) 

(23) 
(24) 
(25) 



(26) 
(27) 
(28) 



Here m| represents the left handed scalar mass squared matricies vplq L and m 2 LL . m 2 R 
represents the right handed scalar mass squared matricies m\j R , rn 2 DR , m? ER and i r n 2 NR . 

3.3 D-term Contributions 



We now consider the D-terms from breaking the Pati-Salam group SU(4) <g) SU(2)^ ® 
SU(2) fl down to the MSSM group SU(3) ® SU(2) L ® U(l)y . These will be family 
independent, but charge dependent, and will pull the six matricies that appear in the 
RGE equations apart. We shall neglect D-term contributions from the broken family 
symmetry which would lead to additional sources of flavour violation. 

The addition to the D-terms have been written down before [201 • The corrections 
are, in matrix notation: 



m 



m z L + glD 2 



QL 



III 



UR 



III 



2 R - {g 2 A - 2g 2 2R )D' 



(29) 
(30) 
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Ill 



DR 



III 



LL 



m 



ER 



in 



NR 



m 



m 2 R -(g 2 + 2g 2 R )D 2 
m 2 L - 3g\D 2 
m R +(3g 2 -2g 2 R )D 2 
m\ + {3g\ + 2g 2 2R )D 2 



in 



h 2 



in 



(31) 
(32) 
(33) 
(34) 
(35) 
(36) 



where in the appendix of Ref. [20], an expression for D 2 in terms of the soft paramaters 
m 2 H and m 2 ^ is derived, 



D 



in 



H 



m 



H 



4\ 2 +2g 2 R + 3gl (37) 
The gauge couplings and mass parameters in EqJSTl are predicted from the model. 
The only free parameter is the coupling is a dimensionless coupling constant which 
enters the potential [20] and should be perturbative. Furthermore, we see that the 
largest that D 2 can be is when is zero, so not only is the order of magnitude of 
D 2 predicted in this model, but we also have an exact upper bound on the value. 



3.4 Soft Gaugino Masses 

The soft gaugino masses are the same as in , which we quote here for completeness. 
The results follow from Eq. (jH2D applied to the SU(4) <g> SU(2) L <g> SU(2) fl gauginos, 
which then mix into the SU(3) ® SU(2) L ® U(l)y gauginos whose masses are given 



by 



M 3 

M 2 
Mi 



\f3mz/2 



T\ + T x ) + [T 2 + T 2 
\Z3m 3 / 2 X Tl 

V3m 3/2 



(T 1 + T 1 )X Tl + (T 2 + T 2 )X T2 



^(T 1 + T 1 )X Tl + ^(T 2 + T 2 )X T2 



(38) 

(39) 
(40) 



l(T 1 + T 1 ) + l(T2 + T 2 ) 
The values of T 4 + Ti and T 2 + T 2 are proportional to the brane gauge couplings g 5l 

and (?5 2 , which are related in a simple way to the MSSM couplings at the unification 
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scale. This is discussed in 

When we run the MSSM gauge couplings up and solve for g 5l and g$ 2 we find that 
approximate gauge coupling unification is achieved by Ti + T\ ^> T 2 + T 2 - Then we 
find the simple approximate result 

Mi « M 3 « M 2 = V3m 3 / 2 X Tl . (41) 
3.5 Soft Trilinear Masses 

So far the soft masses are as in [Hj, with the FN fields and heavy Higgs contribu- 
tions being completely negligible due to the smallness of their F-terms. However 
for the soft trilinear masses these contributions are of order 0(7713/2) despite having 
small F-terms, so FN and Higgs contributions will give very important additional 
contributions beyond those considered in |14j . 

From Appendix[X]we see that the canonically normalised equation for the trilinear 

is: 

A abc = Fj [K T - dj In (k a K b K c ) ] + F m d m In Y abc (42) 

This general form for the trilinear accounts for contributions from non-moduli F- 
terms. These contributions are in general expected to be of the same magnitude as the 
moduli contributions despite the fact that the non-moduli F-terms are much smaller 
|2T] . Specifically, if the Yukawa hierarchy is taken to be generated by a FN field, 9 such 
that ~ Q Vii , then we expect Fg ~ m^^O, and then AA+j = FgdglnYij ~ Vij m z/2 
and so even though these fields are expected to have heavily sub-dominant F-terms 
5 they contribute to the trilinears at the same order 0(m 3 / 2 ) as the moduli, but in a 
flavour off-diagonal way. 

In the specific D-brane model of interest here the general results for soft trilinear 



5 In our model the FN and heavy Higgs vevs are of order the unification scale, compared to the 
moduli vevs which arc of order the Planck scale. 
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masses, including the contributions for general effective Yukawa couplings are given 
in Appendix El From Eqs J4|5l we can read off the effective Yukawa couplings, 

Y hFT hFF = (c)i;(rX z w H 7Z H^ h y a F aa F Px . (43) 

S v ' 

V —P x 

hFFay 

Note the extra group indicies that the effective Yukawa coupling Y hF jT^ has, and 
proper care must be taken of the tensor structure when deriving trilinears from a 
given operator. For Model I and II defined earlier, we can write down the trilinear 
soft masses, A, by substituting the operators in EqElinto the results in Appendix iBl 
Having done this we find the result: 



A = V3m 3 / 7 



d\ + dn + 5dg d\ + dn + 3de d 2 + dn + de 
di + du + 4rfe d\ + dn + 2dg d 2 + dn 
d 3 + dn + 4:dg d 3 + dn + 2de c? 4 



(44) 



where 



d x = X s -X Tl -X T2 (45) 

d 2 = \x s -X Tx - X -X T2 (46) 

d 3 = l -X s -X Tl -X T2 + ]^X T (47) 

d 4 = -X Tl (48) 

d H = (S + S)*X H + (T 3 + T 3 )*X W (49) 

d e = (S + ~B)*(T 3 + T 3 )$X e (50) 



3.6 Why we expect large flavour violation 



According to ref. [22], if the trilinears can be written in a certain manner, then flavour 
violation is expected to be small. However, it is not possible to write the trilinears in 
this manner if either Xn ^ or X^j ^ or Xe ^ 0. It is possible in general to write: 

A ij = Af + Af + S ij (51) 
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If 5ij = 0, then the trilinears factorise. 



.4 



A? 



A* 



A? 



A\ 



Ah 



Ah 



(52) 



If this is true at the SUSY breaking scale, then the FCNC effects are small, and 
the leading order contributions are proportional to Af — Af or Af — Af . If there 
are any contributions which are universal, A , then we can add them in any linear 
combination to Af and Af provided that Af + Af = Af + Af + A . From the 
SUGRA formula for A t ^ it is clear which terms contribute to the universal part, the 
left part and the right part: 



Aij = A hij = F a d a (K - In K h h ) - F a d a In K\ - F a d a In K] +6 V 



(53) 



A R 



We see that 5ij is the terms due to the derivative of the Yukawa. If this is either zero 
or universal, then the A matrix can be written in the restricted form. 

Unfortunately, neither the Higgs contribution or the Froggatt-Nielsen contribution 
can be written in this form. The Higgs contribution to Eql44l is: 



A H oc 



a a a 
a a a 
a a b 



a + {b — a) 





1 



(54) 



The Froggatt-Nielsen contribution is maximally non-universal, and the elements have 
Aij?fc Aki for i k ; j I. In this case we expect there to be the largest contribution 
to flavour violation, assuming that it is not tuned down. ( We could do this either 
by selecting a very small value for the F-term vev by setting Xg m 0, or by setting 
the operator texture in the Yukawas to have very small off-diagonal elements, as the 
A-contribution multiplies Yukawa elements ). Hence we see that the new sources of 
flavour violation will not only contribute to the trilinear terms on at least an equal 
footing as the moduli, but that also they cannot be written in a form where the 
contribtuion to flavour violation is expected to be small. 
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4 Results 



4.1 Benchmark points 

Since the parameter space for this model is reasonably expansive, and the intention 
is to compare different sources of LFV, it is convenient to consider four benchmark 
points, as follows. It should be noted that for all these points, we have taken all 
to be the same, X^ = X?, and also Xh = Xjj. 



Point 


X s 


Xj' 


X H 


Xh 


X e 


A 


0.500 


0.500 


0.000 


0.000 


0.000 


B 


0.536 


0.488 


0.000 


0.000 


0.000 


C 


0.270 


0.270 


0.000 


0.000 


0.841 


D 


0.270 


0.270 


0.578 


0.578 


0.000 



Table 4: The four benchmark points, A-D 



• Point A is referred to as "minimum flavour violation" . At the point Xs = Xj> t 
the scalar mass matrices m 2 are proportional to the identity, and the trilinears 
A are aligned with the Yukawas. Also, if we look back to eq. (J37|) . eq. (J24j) and 
eq. ()25|) we see that for Xs = Xt, which is the case for point A ( and point 
C, and point D ) we see that the upper limit on the magnitude of the D-term 
contribution is zero. As such both m 2 and A will be diagonal in the SCKM 
basis in the abscence of the RH neutrino field. 

• Point B is referred to as "SUGRA". With Xs ^ X^ it represents typical flavour 
violation from the moduli fields; this is the amount of flavour violation that 
would traditionally have been expected with no contribution from the Fh or Fg 
fields. 

• Point C is referred to as "FN flavour violation" . It represents flavour violation 
from the Froggatt-Nielsen sector by itself, without any contribution to flavour 
violation from traditional SUGRA effects since Xs = Xx t as in point A. 
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• Point D is referred to as "Heavy Higgs flavour violation" . It represents flavour 
violation from the heavy Higgs sector, without any contribution from either 
traditional SUGRA effects since X$ = X^, or from FN fields sine Fg = 0. As 
will become apparent, at this point the seesaw mechanism actually helps reduce 
the LFV for fi — * in model I and r — > in both models. 

4.2 Numerical Results 

From the benchmark points defined in Table EJ the F-term vevs were determined, 
and from these the soft paramaters at the high energy scale Mx = 3.10 16 GeV were 
calculated. The soft paramaters were then run down using the 1-loop RGEs of the 
MSSM + v c model. For our numerical results we use a modified version of SOFTSUSY 
. The modifications were made to add the effect of the right-handed neutrino field 
to the RGEs and to decouple them in a manner that allows the neutrino masses and 
mixing angles to be calculated at the low energy scale. As a result of the RGEs having 
to be recoded, all of them are to one loop only in the version that was used here. 

Flavour violation is proportional to non-zero off diagonal elements in the scalar 
mass squared matricies m 2 in the SCKM basis and to non-zero off diagonal elements in 
the trilinears A in the SCKM basis. 6 . Hence, there are two ways to generate flavour 
violation. The first is to have non-zero off-diagonal elements in m 2 of the scalars 
and A at the unification scale. The second is to have non-zero off diagonal elements 
radiatively generated by the /5-function running down to the electroweak scale. It 
is possible to remove the second source by removing the RH neutrino field from the 
model; this allows a disentangling of the see-saw mechanism from the particular source 
of interest, but is unphysical since we know that the neutrinos have to be massive. 

Figure [T] shows numerical results for BR(/i — > ej) for Model I, plotted against the 

gravitino mass m 3 / 2 , where each of the four panels corresponds to each of the four 
6 The SCKM basis is the basis where the yukawas are diagonal at the electroweak scale 
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Figure 1: BR(/i — > erf) for points A-D in model I ( low fi — > e^). The solid line 
represents model I. The dashed line represents an unphysical model with no right- 
handed neutrino field whose purpose is only comparison. The horizontal line is the 
2002 experimental limit from ref. [21] 
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Figure 2: BR(r — > fi^y) for points A-D. The lines coincide in both model I and model 
II. The solid line represents models I and II ( which predict very similar rates for 
r — > fi'y). The dashed line represents an unphysical model with no right-handed 
neutrino field whose purpose is only comparison. The horizontal line is the 2002 
experimental limit from ref. [21] 
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Figure 3: BR(/i — > e^y) for points A-D in model II (high \x — > ej). Note that 
because the rate is much higher the scale is different to that in fig. ^ The solid line 
represents model II. The dashed line represents an unphysical model with no right- 
handed neutrino field whose purpose is only comparison. The horizontal line is the 
2002 experimental limit from ref. [21] 
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benchmark points A-D. As 7713/2 increases the sparticle spectrum becomes heavier. 
This will look different at each paramater point, but the physical masses are expected 
to be of the same order of magnitude as the gravitino mass. As such, high gravitino 
masses will start to reintroduce the fine-tuning problem of the gluino mass being too 
high. Point A corresponds to minimum flavour violation, where the only source of 
LFV is from the see-saw mechanism, which for Model I is well below the experimental 
limit, shown as the faint horizontal dashed line. Point B has LFV arising from 
SUGRA, with the FN and heavy Higgs sources of LFV switched off, and in this case 
we also show the results with the see-saw mechanism switched off (dashed curve) as 
well as with the see-saw mechanism with SUGRA contributions to LFV (solid curve). 
In both cases the results are below the experimental limit for 7773/2 above 500 GeV. 
Point C is the FN benchmark point, and for this case we see that the experimental 
limit is violated over the entire range of 777,3/2 shown, with the see-saw mechanism 
making very little difference. Point D shows the heavy Higgs point, for which the 
experimental limit is violated for 7773/2 below 1000 GeV. Interestingly, the effect of 
switching off the see-saw mechanism in this case (dashed curve) is to increase the rate 
for BR(/i — > ej). 

Figure El shows results for BR(r — > fi'y) for Model I, plotted against the gravitino 
mass 777,3/2. Point A for minimum flavour violation is below the experimental limit, 
as is point B corresponding to SUGRA, with the see-saw mechanism switched off 
corresponding as before to the dashed curve. Point C corresponding to FN violates 
the experimental limit for lower 7773/2, with a rather large effect coming from the see- 
saw mechanism. Point D shows the heavy Higgs point, with the effect of the see-saw 
mechanism being to reduce BR(r — > /ry) in conjunction with the LFV coming from 
heavy Higgs, similar to the analagous effect observed previously. 

Figure El shows the analagous results for BR(yU — > ej) for Model II. As expected 
model II, which is supposed to give a high rate for ji — > ej, does give results close 
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to the experimental limit for points A and B, and the limit is now well exceeded for 
points C and D. By increasing the gravitino mass sufficiently (which increases all 
the sparticle masses) it is possible to respect the current experimental limit, but at 
the expense of a very heavy superpartner spectrum. However it is worth noting that 
for benchmark points C,D the value of Xt is almost half its value corresponding to 
points A,B. According to EqJ41lthis implies that for points C,D the gaugino masses are 
almost half their values corresponding to points A,B, leading to reduced fine-tuning 
for a given 7713/2- 

The results for BR(r — ► fij) for Model II are almost identical to those shown for 
Model I in Figure El which is as expected since the only difference between the two 
models is in the 13 element of the neutrino Yukawa matrix. 

5 Conclusions 

We have investigated the phenomenological impact of different sources of lepton 
flavour violation arising from realistic D-brane inspired models based on supergravity 
mediated supersymmetry breaking, where the origin of flavour is due to Froggatt- 
Nielsen (FN) operators. We have discussed four distinct sources of lepton flavour 
violation in such models: minimum flavour violation, arising from neutrino masses 
and the see-saw mechanism with renormalisation group (RG) running; supergravity 
flavour violation due to the non-universal structure of the supergravity model; FN 
flavour violation due to the FN fields developing supersymmetry breaking F-terms 
and contributing in an non-universal way to soft trilinear terms; and finally heavy 
Higgs flavour violation arising from the Higgs fields used to break the unified gauge 
symmetry which may behave analagously to the FN fields. 

In order to quantify the relative effects, we studied a particular type I string 
inspired model based on a supersymmetric Pati-Salam model arising from intersecting 
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D-branes as proposed in [T3], but here supplemented by a £7(1) family symmetry 
with the quarks and leptons described by the set of FN operators as in [12] • We 
have derived the soft supersymmetry breaking masses for the model, including the 
new flavour vilating contributions to the soft trilinear masses arising from the FN 
and heavy Higgs fields. We then performed a numerical analysis of LFV for four 
benchmark points, each chosen to highlight a particular source of flavour violation, 
with the benchmark points C and D coresponding to LFV arising from the FN and 
heavy Higgs fields giving by far the largest effects. Since the new contributions are 
dominantly from the trilinears A, the amount of flavour violation is therefore strongly 
dependent on the choice of Yukawa matrices at the unification scale. For example 
the huge difference in the rate of \x — > between Model I and Model II is simply 
generated by changing the (1,3) element of Y u . Also we find that \x — > e'j is more 
constraining than r — > fj,j. 

The numerical results show that the contributions to LFV from Yukawa operators 
with the heavy Higgs sector and the Froggatt-Nielsen sector can give the dominant 
contributions to LFV processes, greatly exceeding contributions from SUGRA and the 
see-saw mechanism, and should be taken into account when performing phenomeno- 
logical analyses of supergravity models. 
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A Soft terms from supergravity 

We summarise here the standard way of getting soft SUSY breaking terms from 
supergravity. Supergravity is defined in terms of a Kahler function, G, of chiral 
superfields (0 = h,C a ). Taking the view that the supergravity is the low energy 
effective field theory limit of a string theory, the hidden sector fields h are taken to 
correspond to closed string moduli states (h = S,Ti), and the matter states C a are 
taken to correspond to open string states. In string theory, the ends of the open 
string states are believed to be constrained to lie on extended solitonic objects called 
Dp— branes. 

Using natrual units: 



K(<f>, 4>) is the Kahler potential, a real function of chiral superfields. This may be 
expaned in powers of C a : 



Kab is the Kahler metric. W(4>) is the superpotential, a holomorphic function of 
chiral superfields: 



We expect the supersymmetry to be broken; if it is broken, then the auxilliary 
fields Ftj, 7^ for some 0. Lacking a model of SUSY breaking we can proceed no 
further without paramaterising our ignorance. We do this using goldstino angles. We 

introduce a matrix, P that canonically normalises the Kahler metric, P^K-jjP = 1 7 

7 The subscripts on the Kahler potential Kj means djK. However, the subscripts on the F-terms 
are just labels. 




(55) 



K = K(h, h) + K m {h, h)CaC h + [-Z ab (h, h)C a C b + /i.c.J + ... 



(56) 



W = w(h) + -fi ab (h)C a C b + -Y abc C a C b C c + ... 
z 



(57) 
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. We also introduce a column vector G which satisfies 0^0 = 1. We are completely 
free to paramaterise G in any way which satisfies this constraint. 

Then the un-normalised soft terms and trilinears appear in the soft SUGRA break- 
ing potential [2"o] : 

V soit = ml b CaC b + {^A abc Y abc C a C b C c + h.c)j + ... (58) 
The non-canonically normalised soft trilinears are then: 

A-abc^abc TTtTi ^ K m Y abc -\- d m Yabc \\K ) 3 m Ke a Yd bc 

\W\ 1 vv ' 

+ (a^b) + (a^ c))l (59) 



In this equation, it should be noted that the index m runs over h, C. However, 
by definition, the hidden sector part of the Kahler potential and the Kahler metrics 
are independent of the matter fields. 

Assuming that the terms dcY abc ^ 0, the canonically normalised equation for the 
trilinear is: 

abc 



A abc = Fj [Kj - dj In (K a K b K c ) \ + F m d m In Y abc (60) 

If the Yukawa heirarchy is taken to be generated by a Froggatt-Nielsen field, <f) such 
that Y oc <fr p , then we expect F^ oc 7Ti 3 / 2 0, and then F^d^lnY oc m 3 / 2 and so even 
though these fields are expected to have heavily sub-dominant F-terms, they con- 
tribute to the trilinears on an equal footing as the moduli. 

If the Kahler metric is diagonal and non-canonical, then the canonically normalised 
scalar mass-squareds are given by 



ml^-FjFjdjd! (hi K a ) (61) 



And the gaugino masses are given by 



M » = vkr/ M - (62) 
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Where f a is the 'gauge kinetic function', a enumerates D-branes in the model. In 
type I string models without twisted moduli these have the form fg = S; f& = T\ 

Specifically, we use a Kahler potential that doesn't have any twisted-moduli \27\ : 



K 



In iS + S - Cf 1 

In (t 2 - T 2 - 

|C 5l52 | 2 



ci 2 



In Ti + Tx- C\ 



cl 



[T 3 -T 3 

|2 



c 2 51 



Cf 1 



|C 951 



/ _\ 1/2 / _ 

(S + S) (T3 + T3 

| C 95 2 |2 



+ 



1/2 



1/2/ — \ 1/2 

T 2 + T 2 ) (T3 + T3, 



7^ \V2 



(Tx + Tx) ' (T 3 + T 3 



7^ \V2 



(63) 



The notation is that the field theory scalars, the dilaton S and the untwisted 
moduli Ti originate from closed strings. Open string states C\ are required to have 
their ends localised onto D— branes. The upper index then specifies which brane(s) 
their ends are located on, and if both ends are on the same brane, the lower index 
specifies which pair of compacitified extra dimensions the string is free to vibrate in. 

B Paramaterised trilinears for the 42241 Model 



We here write the general form of the trilinear paramaters assuming nothing 
about the form of the Yukawa matricies. 



A r 

cf 1 C B i 5 aC7 5 i s 2 



V3m 3/2 {jf, 
+ X Tl 
+ X T2 



l + [S + S)d s lnY abc 
l + ^ + Tx^lnlW 
l + (T 2 + T 2 )d T2 \nY abc 



+ X T3 (T 3 + T 3 )d T JnY c 



abc 



X H S + S) 2 Hd H InY 



abc 



+ X-u T3 + T3 ) 2 Hd w \nY abc 



+ X e (S + Sy [T 3 + T 3 ) 1 6d e In Y abc 



(64) 
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A c^c^c^2 = V3m 3/2 lx, 



S + S)d s lnY abc 



+ X Tl -1+ Tj + Ti )d Tl lnY c 



abc 



+ X T JT 2 + T 2 )d T AnY c 



abc 



+ x. 



T 3 



■-{T 3 + T 3 )d T3 \nY abc 



+ X H [S + S) 2 Hd H \nY ahc 

+ X 1T (T 3 + T 3 ) h Hd w \nY abc 

+ X e (S + S) 1 (T 3 + T 3 ) 1 #<9 e In F abc 



\/3m 3 / 2 <^ X 



5 + 5 9 s i n y c 



«6c 



+ x Tl 


-1 4 


-(21- 


hT 




in y afec 


+ x T2 


-1 4 


(r 2 - 


hT 2 ; 


d T2 


In Y abc 



+ x 



r 3 



T 3 + T 3 )a T3 lny c 



abc 



+ X e ( 1 S + 5) 5 (T 3 + T 3 ) 5 ^ e lnr c 



.)^%inr afec 



abc 



^r> 5 i r> 5 i r> 5 i 

°1 °2 °3 



+ X T2 {T 2 + T 2 )d T2 \nY abc 
+ X T3 (T 3 + T 3 )d T3 \nY abc 
+ X H (s + sf Hd H \nY abc 
+ X 7T (T 3 + T 3 ) ll Hd w \nY c 
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abc 



+ X (S + S) 1 (T 3 + T 3 ) 1 Ode In Y abc } (67) 

C n = 1 operators 

The n = 1 Dirac operators are the complete set of all opearators that can be con- 
structed from the quintilinear FFhHH by all possible group theoretical contractions 
of the indicies in 

Off = F aa F Px hlH, z H^ (68) 

We define some SU(4) invariant tensors C and some SU(2) invariant tensors R as 
follows 8 : 



(CO? 


eft 

= Sp 




t a/3u;x 


l°ioJ P7 










= 5 v 







(69) 



Then the six independent SU(4) structures are: 



A. 


(Cif a (Gx)J 




B. 






C. 


^ \wx (n \P r < 


= 8(^*2 -^?) 


D. 




= 2(^ + 5^) 


E. 


{Cif p mi 




F. 







(70) 



8 The subscript denotes the dimension of the representation they can create from multiplying 4 
or 4 with 4 or 4. For example (Ci 5 )^4 7 4' > = 15^ 
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Operator Name 


Operator Name in HSI 


QUh 2 


QDhi 


LEhi 


LNh 2 


U 


U 


1 


-1 
1 


-1 
1 


1 


U 


U 


1 






1 


qAc 


QM 












QAd 


QT 


2V2 
5 

V2 


V2 
5 


V2 
5 


2V2 
5 

v/2 


QAe 


o v 








Af 


o u 


5 


2V2 

5 


2V2 

5 


5 


QBa 


o c 


1 

V5 


1 

V5 


-3 

VE 


-3 


nBb 

u 


(J 


1 


-1 


-3 


3 


VE 


Vf> 






U 


U 


U 


[2 

V 5 


" 3 V5 


U 


QBd 


X 


2V2 
5 


V2 
5 


-3v/2 
5 


-6\/2 
5 


QBe 


o z 


/ ^ 
V 5 








6 y 5 


Bf 


Y 


V 2 
5 

V2 


2v2 

~E~ 

V2 


— 6\/2 
5 


— 3v2 
5 














(J 


n F 

yj 


V l 


—y l 


n 
u 


n 
u 


qCc 


qE 


n 




Zj 


n 


n 


QCd 


b 


4 

VE 


2 

V5 








QCe 


QN 


2 











nCf 




2 


4 

V5 


n 


n 


o 


r^d 
O 


pi 

V 5 


j2 
V 5 


2 vl 


2 vl 


QDb 


o e 


[2 

V 5 


[2 
V 5 


-2,/^ 

Z V 5 




QDc 


0° 





2 

V5 


4 

^5 





QDd 


H 


4 
5 


2 
5 


4 
5 


8 
5 


QDe 


0° 


2 

V5 


n 


n 


4 


Df 


of 


2 
5 


4 
5 


8 
5 


4 

5 

V2 


QEa 


0° 








QEb 


Qh 








-V2 


V2 


O 


s~\ j 

O l 








2 





QEd 


o j 








2 

Vf> 


4 


QEe 


o 1 











2 


Ef 


J 








4 


2 

v/5 


QFa 


o p 


AV2 


AV2 


3^ 


3^ 


5 


5 


5 


5 


QEb 


Q 


AV2 
5 


-4V2 
5 


-3V2 
5 


3V2 
5 


qFc 


qR 





8 
5 


6 
5 





QFd 


L 


16 


8 


6 


12 


5 f 


5^ 




5^ 


QFe 


qK 


5 








6 
5 


QFf 


s 


8 

bVE 


16 

bVE 


12 

5^ 


6 

hVE 



Table 5: Operator names, CGCs and names in ^B] 
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And the six SU(2) structures are: 



w y <2 y w 

d. e ws e°* (R 3 ) s y l (R 3 ) z t r q = 6*J* - l -e wy e™ 

(RiY y (Ri)l = ^51 
f- (^3)^(^3):; = (71) 

All possible n = 1 operators were then named A ...0 z a ...0 j in [IE]. We rename 
them here in a manner consistent with the n > 1 operators O^'-*, so that the names 
are O nn where IT is the SU(4) structure and 7r is the SU2 structure. See Table El for 
the translation into the names of ref.^Hl an d the CGCs. 

All of these operators are operators for the case without a U(l) family symmetry. 
In the case when there is, we follow the prescription: 

"^°"{m- x ) (72 > 

Where pu = \Xojj\ is the modulus of the charge of the operator. If the charge 
of the operator is negative, then the field 9 should be replaced by the field 9. The 
prescription makes the operator chargeless under the U(1)f while simultaneously not 
changing the dimension. 

D n > 1 operators 

In the case that n > 1, there will be more indicies to contract, which allows more 
representations, and hence more Clebsch coefficients. To generalise the notation, 
it is necessary only to construct the new tensors which create the new structures. 
However, it will always be possible to contract the new indicies between the H and H 
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fields to create a singlet HH which has a Clebsch of 1 in each sector u, d, e, v. In this 
case, the first structures are the same as the old structures, but with extra 5 symbols 
which construct the HH singlet. 

Thus taking a. n = 2 operator, say ,Fb , which forms a representation that could 
have been attained by a n = 1 operator, the Clebsch coefficients are the same. This 
is what we mean by (9 n/ri7r , as we have only used n > 1 coefficients which are in the 
subset that have n = 1 analogues. 
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